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Gluon field configurations with nonzero topological charge generate chirality, inducing V- and CP- 
odd effects. When a magnetic field is applied to a system with nonzero chirality, an electromagnetic 
current is generated along the direction of the magnetic field. The induced current is equal to the 
Chiral Magnetic conductivity times the magnetic field. In this article we will compute the Chiral 
Magnetic conductivity of a high-temperature plasma for nonzero frequencies. This allows us to 
discuss the effects of time-dependent magnetic fields, such as produced in heavy ion collisions, on 
chirally asymmetric systems. 



I. INTRODUCTION 

Quantum chromodynamics (QCD) contains gauge field 
configurations which carry topological charge [l| . These 
configurations interpolate between the different vacua of 
the gluonic sector of QCD Q and induce interesting V- 
and CV-odd effects 0]. Neglecting the masses of quarks 
(as appropriate at high energies), it holds for each indi- 
vidual flavor that 4] 



AiV5 = A{Nr - Nl) 



(1) 



where AN5 denotes the change in chirality {N5) which is 
the difference between the number of modes with right- 
and left-handed chirality. In the limit of zero quark mass 
A^5 is also equal to the total number of particles plus 
antiparticles with right-handed helicity minus the total 
number of particles plus antiparticles with left-handed 
helicity. Right-handed helicity means that spin and mo- 
mentum are parallel, whereas left-handed helicity implies 
they are opposite. 

There exist different mechanisms to generate topolog- 
ical charge during heavy ion collisions. One possible way 
is by longitudinal fields created just after the collision 
[liSiHi another due to QCD sjAaleron transitions in the 
quark-gluon-plasma [1, (9|, [l^, [O] , and also plasma insta- 
bilities can lead to generation of topological charge p^ . 
Furthermore, it has been argued that instanton "ladders" 
may describe a significant fraction of multi-particle pro- 
duction at high energies [H, [ij, [l^ . Finally, metastable 
V— and CV— odd domains may exist in the quark-gluon 
plasma close to the critical temperature flG^. No mat- 
ter what the precise mechanism behind the generation of 
topological charge is, this leads to generation of chirality 
as can be seen from Eq. ([1]). 

Assuming that the 9 angle vanishes there is no ex- 
plicit V- and CT'-breaking in QCD. Hence positive and 
negative topological charges are being generated with 
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equal probability. But because of fluctuations a finite 
amount of topological charge can be generated in each 
event within a given region of phase space determined by 
the experimental acceptance. Only if one averages over 
many events the net produced topological charge should 
vanish. Therefore in each event a difference between the 
total number of right- and left-handed particles can be 
expected. 

When two heavy ions collide with nonzero impact pa- 
rameter, a (electromagnetic) magnetic field of enormous 
magnitude is created in the direction of angular momen- 
tum of the collision; it has been evaluated in Refs [13, [ill 
(see also [l^ for a proposal to utilize electromagnetic 
fields in the search for disoriented chiral condensates) . If 
a nonzero chirality is present in such a situation, an elec- 
tromagnetic current will be induced in the direction of the 
magnetic field. This is the so-called Chiral Magnetic Ef- 
fect [H HO, HI HI; see recent Ref. [H] for the first study 
of this effect in lattice QCD. To understand this effect 
qualitatively let us imagine a situation in which the V- 
and CP-odd processes made positive, so that we have 
an excess of quarks plus antiquarks with right-handed 
helicity. In a background magnetic field, the quarks will 
align their magnetic moments along the magnetic field. 
And assuming the quarks can be treated as massless, the 
momenta of the quarks will be aligned along the field 
as well. Consequently a quark with left-handed helicity 
tends to move exactly in the opposite direction to a quark 
with right-handed helicity. Since the magnetic moment 
is proportional to the charge, an antiquark with right- 
handed helicity will move exactly opposite to a quark 
with right-handed helicity. Accordingly, in this case of 
positive an excess of positive charge will move paral- 
lel to the magnetic field and an excess of negative charge 
will go in the opposite direction. Thus an electromag- 
netic current is generated along the magnetic field. 

In a heavy ion collision this current leads to an excess 
of positive charge on one side of the reaction plane (the 
plane in which the beam axis and the impact parame- 
ter lies) and negative charge on the other; the resulting 
charge asymmetry is also modulated by the radial fiow 
and the transport properties of the medium. This charge 
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asymmetry can be investigated experimentally using the 
observables proposed in Ref. 0| . Preliminary data of the 
STAR collaboration has been presented in Refs. [25112^. 
Implications of the Chiral Magnetic Effect on astrophys- 
ical phenomena have recently been discussed in Ref. ,27il ; 
another astrophysical implication can be found in [28| . 

A system of massless fermions with nonzero chirality 
can be described by a chiral chemical potential ^5 which 
couples to the zero component of the axial vector current 
in the Lagrangian. The induced current in such situa- 
tion can be written as j — a^B, where cr^ is the Chiral 
Magnetic conductivity. For constant and homogeneous 
magnetic fields its value is determined by the electro- 
magnetic axial anomaly and for one flavor and one color 
equal to [H, HI [131 (see also [IH) 



and quark gluon plasma are equal (up to a trivial factor of 
NcJ2fl})- Since we do not take into account higher or- 
der corrections, some of our results for QCD will only be 
a good approximation in the limit of very high tempera- 
tures where the strong coupling constant as is sufficiently 
small. 

We will take the metric g'^'^ = diag(-|-, — , — , — ). 
The gamma matrices in the complete article satisfy 
{7*^1 7"} = 25'^". We will use the notation p for both 
the four- vector = {po,p) and the length of a three- 
vector p = \p\. 



II. KUBO FORMULA FOR CHIRAL 
MAGNETIC CONDUCTIVITY 



a^{uj = 0,p = 0) = ao 



27r2 



Ms, 



(2) 



where u and p denote frequency and momentum respec- 
tively, and e equals the unit charge. For a finite number 
of colors Nc and flavors / one has to multiply this result 
by ^cJ2f If where qj denotes the charge of a quark in 
units of e. The generation of currents due to the anomaly 
in background fields or rota ting sys tems is also discussed 
in related contexts in Refs. [U [H, [H, [H, [111 . 

For constant magnetic fields which are inhomogeneous 
in the plane transverse to the field one finds that the total 
current J along B equals ^22j , 



J : 



27rJ TT ' 



(3) 



where is the length of the system in the z-direction 
and the flux $ is equal to the integral of the magnetic 
field over the transverse plane. 



(4) 



The floor function lx\ is the largest integer smaller than 
X. The quantity [e$/(27r)J in Eq. ^ is equal to the 
number of zero modes in the magnetic field !30|. 

To compute the current generated by a configuration 
of specific topological charge, one should express fi5 in 
terms of the chirality . By using the anomaly relation 
one can then relate to the topological charge. This is 
discussed in detail in Ref. [2^ . 

The aim of this paper is to study how a system with 
constant nonzero chirality responds to a time dependent 
magnetic field. This is interesting for phenomenology 
since the magnetic field produced with heavy ion colli- 
sions depends strongly on time. To obtain the induced 
current in a time-dependent magnetic field, we will com- 
pute the Chiral Magnetic conductivity for nonzero fre- 
quencies and nonzero momenta using linear response the- 
ory. We will compute the leading order conductivity and 
leave the inclusion of corrections due to photon and or 
gluon exchange for future work. In leading order the Chi- 
ral Magnetic conductivity for an electromagnetic plasma 



For small magnetic fields, the induced vector current 
can be found using the Kubo formula. This formula tells 
us that to first order in the time-dependent perturbation, 
the induced vector current is equal to retarded correlator 
of the vector current with the perturbation evaluated in 
equilibrium. More explicitly, one finds that 



d^x' Il^''{x,x')A,{x'), 



(5) 



where j^{x) — e'ip{x)^f^ip{x) and the retarded response 
function 11^^^ is given by 

U^''{x,x')=t{[r{x),,f{x')])9{t~t'). (6) 

The equilibrium Hamiltonian is invariant under trans- 
lations in time and space, therefore we can use that 
Il'^'^{x, x') — li'^ [x — x'). Let us take a vector field of the 
following specific form A^{x) = A^{p)e^^P^ . The Kubo 
formula now becomes, 



where 



{rix))^Il^^ip)AAp)e-^P'', 



il'^'ip)^ d^xe^P^U^^^ix). 



(7) 
(8) 



In order to compute the Chiral Magnetic conductiv- 
ity we will take a time-dependent magnetic field pointing 
in the z-direction. Because of Faraday's law (V x E = 
—dB/dt), such time-dependent magnetic field comes al- 
ways together with a perpendicular electric field. Let 
us choose a gauge such that the only component of the 
vector field that is non- vanishing is Ay. Then Bz{x) = 
dxAy{x) so that Bz{p) = ip^A^{p). Using Eq. jT]) we 
find that the induced vector current in the direction of 
the magnetic field can now be written as 



(J,(x))=a^(p)B,(p)e-^^ 
where the Chiral Magnetic conductivity equals 



1 ~ 



ip^ 



2ip'- ^ ^ ' 



(9) 



(10) 
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The right-hand side of the last equation is a result of 
rotational and gauge invariance. We can now write the 
Chiral Magnetic conductivity in the following way 



with 



(11) 



(12) 



The last two equations show that the Chiral Mag- 
netic conductivity follows directly from evaluating (p) , 
which is the spatially antisymmetric part of the off- 
diagonal retarded current-current correlator, or equiva- 
lently the photon polarization tensor. 

Another related quantity that follows from the off- 
diagonal part of the photon polarization tensor is the Hall 
conductivity. The Hall current is generated in the pres- 
ence of a magnetic field that is perpendicular to an elec- 
tric field. Unlike the Chiral Magnetic current, the Hall 
current is in a direction perpendicular to the magnetic 
field. In order to obtain the Hall conductivity one usually 
computes the photon polarization tensor in the presence 
of a homogeneous background field. The electric field is 
then treated as a perturbation (see e.g. Ref. [37.] for a 
recent calculation using a holographic model of QCD). 

In general, the Chiral Magnetic conductivity will be 
complex. Let us therefore write 

^x(p) = <^x(P) + i^x^P)^ (13) 
where cr^(p) and o'^(p) are real functions given by 



a'^ip) = -lmG^(p), 
<(P) = -^ReGj,(p). 



(14) 
(15) 



For convenience we will write the zero momentum limit 
of the Chiral Magnetic conductivity as follows 



: UJ,p). 



(16) 



The real part of the conductivity is an even function of lj 
while the imaginary part is odd, i.e. o'^('-^) = ^'xi~^) ^^"^ 
CT^'(cj) = — ct"(— w). If we apply a homogeneous magnetic 
field B{t) — cos{Lot)z the response j{t) = {jz{t)) will 
be 

j{t) = [<7'^{lj) cos{u;t) + a'^{ij) sm{ujt)] (17) 

Hence the real part and imaginary part correspond to 
the in- and out-phase response respectively. To find the 
response to a general time-dependent magnetic field B = 
B{t)z one first has to compute the Fourier transform of 
the magnetic field 



B{lu) = / dte'^^^Bit). 



(18) 



The response will then be 



j{t) = — [a'[uj) cos(wt) + ct"(cj) sin(wt)] B{uj). 
Jq tt 

(19) 

The real and imaginary parts of the Chiral Mag- 
netic conductivity related to each other by the Kramers- 
Kroning relation 



1 



dqo 



-oo 90 - w 

1 f°° 
a'^{u) = V dgo 

^ J —OO 



(20) 



(21) 



qo-u} 

In the next section we will study the ferniion prop- 
agator in the presence of a chiral chemical potential. 
Then using this propagator we will compute the retarded 
current-current correlator from which we can derive the 
Chiral Magnetic conductivity. 



III. FERMION PROPAGATOR 

The bare fermion propagator as a function of Euclidean 
momentum Q in the presence of a chiral chemical poten- 
tial equals 

S{Q) - -^jp , (22) 

here (D„j — (2m + 1)ttT with to G Z is a fermionic Mat- 
subara mode. By computing the inverse, this propagator 
can be written as 



(23) 



where we have defined Qj. = (iijJ± , q) with ili± = ujrn — 
and /i± = fizL^^. Accordingly Q±^ = {iuj±, —q) and 
Q\ = -{Cj'i^ + q^). Furthermore ^± = Q±"f^. The right- 
{+) and left-handed (— ) chirality projection operators 
are given hyV± = ^(1 ±75). They satisfy V"^ = V± and 
V+V- = 0. Let us now define 



A±((7o,g) 



qoTEq' 



(24) 



where Eg = |q|. Furthermore we introduce q = q/\q\ 
and (jj. = (l,±q). With these definitions the fermion 
propagator can be rewritten as 



5(Q) = i M^^s,q)VsJ^q^. 

s,t=± 



(25) 



As can be seen from Eq. the propagator consists of 
two parts describing the modes with right- (s — +) and 
left-handed chirality (s = — ). Because we took m = 
the opposite chiralities do not mix. The different values 
of t correspond to particles (modes with positive energy, 
t = +) and antiparticles (modes with negative energy, t = 
— ). In the following section we will use this expression for 
the propagator to compute the retarded current-current 
correlator. 
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IV. COMPUTATION OF RETARDED 
CORRELATOR 



which amounts to replacing iujn by po + is in Eq. (I3ip 
We obtain 



We will compute the retarded correlator de- 
fined in Eq. (jl2p . using the imaginary time formalism 
of thermal field theory. The retarded correlator can be 
obtained from the Euclidean correlator by analytic con- 
tinuation in the following way 



(26) 



where e = 0+. 

At very high temperatures, the gluons and quark 
masses can be neglected to first approximation and the 
current-current correlator is a convolution of two bare 
massless fermion propagators S{Q) (see e.g. Ref. [38,] for 
the correct expression). Using Eq. (|12p we find 



e 

2(3 



E 



d3g 

(2^ 



e*-'''=tr [7^5(Q)7-''5(P-t-Q)] . 

(27) 

With use of Eq. ([25]) and the properties of the chirality 
projection operators the integrand of Eq. (|27[) can now 
be written as 
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At(ia)s,q)A„(ia;s + ii^n.p + q)q'^{p + q)"^. (28) 

It can be seen from the last equation that the opposite 
chiralities (s — ±) do not mix. As long as m = 0, the 
Chiral Magnetic conductivity is a sum of a contribution 
from purely right-handed modes and purely left-handed 
modes. 

We can now use that tr[7''7''7''7'^7^] = -Aie^'•'P'' 
where e^'''"^ is the complete antisymmetric tensor with 
g0i23 ^ rj,^^^ follows that 



As a result we obtain 

.2_ r 



3i {a'}P - a%') . (29) 



GUP) 



E 



(27r)3 



E ^ 

s.i.u— ± 



-I- 



'P+9 



At(ia)s, q)A„(ia)s -I- iw„,p -I- q). (30) 



We can now perform the sum over Matsubara frequen- 
cies. Using that cj„ = 2miT is a bosonic Matsubara 
frequency one finds 

i^At(ja;s,q)A„(ia;5 +ibJn,P + q) = 



tfl{Eq - t^s) - Un{Ep+q - Ufls) + ^{u - t) 



tEo - uE, 



(31) 



■p+9 



where h{x) = [exp(/3a;)-f 1]^^ is the Fermi-Dirac distribu- 
tion function. We can now perform the analytic contin- 
uation in order to obtain the retarded correlator (p) , 



(27r)3 E • 

^ ' s,t,u=± 



9' p' 

t u- 

E„ 



E.„ 



tn{Eq ~ tfls) - Un{Ep+q - U^,s) + i(li - t) 



Pa 



tEq — uEp+q 



(32) 



Before we compute the integral over q let us try to inter- 
pret Eq. ([32]) . The retarded correlator G|^(p) has a real 
part when pq — uEp^q — tEq. From Eq. ([TS]) it can be 
seen that in that case the Chiral Magnetic conductivity 
acquires an imaginary part. When po — uEp^q — tEq 
the virtual particles in the loop of the photon polariza- 
tion tensor go on the mass-shell, which by the optical 
theorem corresponds to production or scattering of real 
particles in the electromagnetic field. 

If u = — t this implies that particle antiparticle pairs 
are produced from the electromagnetic field that oscil- 
lates with frequency po- If t ~ u particles {u — t — 1) 
or antiparticles [u = t — ~1) scatter from the electro- 
magnetic field and acquire or lose some momentum and 
energy. Let us take a closer look at the pair production 
process. 

If Po > the produced particles and antiparticles have 
energies Epj^q {u = 1) and Eq {t = —1) respectively. If 
the system we consider consists mainly of particles, fis 
is positive. At zero temperature all particle states up to 
the Fermi energy fj,s are filled, so then it is impossible 
to produce particles with energy per particle less than 
fis due to Pauli blocking. This is reflected in the Fermi- 
Dirac distributions in the integrand of Eq. (|32l) , for u = 
t — —1, [Is > 0, and r = 0, the integrand vanishes if 
Ep^q < fis- If Ep-^-q > fis there is no Pauli blocking, in 
that case Eq > — p. Hence at T = 0, G^(p) develops 
a real part and particle-antiparticle pairs with chirality 
s can be produced if pq > 2/is — p. If both Ep^q and Eq 
are larger than fig the produced pair does not feel the 
influence of the nonzero chirality, so then the imaginary 
part of the Chiral Magnetic conductivity should vanish 
as well. In that case Eq < so that Ep+q < fj.^ + p 
which gives po < 2/1^ -I- p. So concluding we expect the 
Chiral Magnetic conductivity to have an imaginary part 
for 2fis ^ p < Pq < 2/is +P due to pair production. 

Now let us continue the calculation of G|{(p). The 
last term in the integrand of Eq. ([5^ . vanishes after 
summing over chiralities (s). It then follows that the 
function G}^{p) is ultraviolet finite because the high- 
momentum part of the integrand is exponentially sup- 
pressed by the Fermi-Dirac distribution functions. We 
can therefore shift and refiect the integration variable as 
follows q — > —q — p in order to make both Fermi-Dirac 
distribution functions dependent on Eq. Then after in- 
terchanging t with u and inversions t —t and u — s- —u 
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we arrive at 



le 
~2 



(2^ 



p+q 



fl{Eq - fis) - n{Eq + /i^) 

Po + + tEq + wi?p+g 
Here we used that since i = ±, 



(33) 



t [n{Eq -tfis) - n{Eq + t^s)] = 

+ ^,). (34) 

We can now perform the sum over u which gives us 



(27r) 



E 



s,t=± 

n{Eq - /is) - n(£'q + /i^) 



(35) 



(po + + t-Bg)2 - El+q 

Now it is possible to perform the angular integral. We 
use that the fact that the integral has to be proportional 
to p\ Hence we can replace by [q ■ p)p^ /p^ ■ We find 
after integrating over angles 



- 2 i 7 

le p p 



pI 



d(?/((?)x 



{2q + tpo) log 

t=± 
where 



IGtt- p Jo 

{pQ + ie + tqf - {q+pf 
{po + ie + tqY - (q - pY 



/(?) = E * ^ ^'^^ ^ + 



(36) 



(37) 



s=± 



We now have all results in order to obtain the Chi- 
ral Magnetic conductivity which follows from combining 
Eq. ((Ml) with Eq. (HS]). We wiU discuss the result in the 
next section. 



V. COMPUTATION OF CHIRAL MAGNETIC 
CONDUCTIVITY 

We will discuss the leading order contribution to the 
Chiral Magnetic conductivity as a function of frequency 
and momentum. The zero frequency and momentum 
value is constrained by the axial anomaly, so loop correc- 
tions by gluons and/or photons will not alter this result 
[121 • However, loop corrections will change the conduc- 
tivity at nonzero frequencies. 

In an electromagnetic plasma, the leading order result 
is a sensible approximation since loop effects are of order 
q;em which is small. The results for an electromagnetic 
plasma are only reliable if the temperature is larger than 
the mass of the electron, since we have assumed massless 
particles. In a quark gluon plasma loop corrections are 
only negligible at very high temperatures due to asymp- 
totic freedom. Hence one should keep in mind that the 
leading order result is in QCD only a valid approximation 
at high temperatures. 



A. Zero frequency limit 

Let us first rederive the zero frequency, zero momen- 
tum limit of the Chiral Magnetic conductivity. Since 



lim lim y log 



t=± 



{po + i€ + tqf - {q + pf 
{Pq + ie + tqY - {q- pY 



2p 

q 



it immediately follows from Eq. (j36|l that 



limGj,(po = 0,p) = ^p' 



(38) 



(39) 



As a result we recover the known zero frequency result of 
the Chiral Magnetic conductivity for homogeneous mag- 
netic fields (p — ^ 0), 



CTo = a^{uj = 0) = ^M5- 



(40) 



The zero frequency limit is independent of fi and T. Since 
this value is constrained by the axial anomaly, loop cor- 
rections by gluons and/or photons will not alter this re- 
sult fH. 



B. Imaginary part 

To obtain the imaginary part of the Chiral Magnetic 
conductivity we need to compute the imaginary part of 
the logarithm in Eq. (j36p . This is a sum of step functions 
times TT. We find using that q > and p = |p| > 0, 



Im^(2g-f ipo)log 



t=± 



{pQ + ie + tqY -{q + pf 
(po + ie + tqY - (g - pf 



IT [2q - \po\e{pl ~ p2)] [e{q+ -q)- 9{q- - q)] 

+ irpoOip^ - pI) [eiq - q+) + e{q - g_)] , (41) 

where q± = i|po ±p|. 

To compute the imaginary part of the Chiral Magnetic 
conductivity we furthermore need the to perform the fol- 
lowing two indefinite integrals 



dqf{q)=T St log [l + e('+*^»)/^ 



s,t=± 



I dqqfiq) = ^ st{qTlog 
•' s,t=± 



-T^Lia 



^Jq+tti.s)/T 



}, (42) 



here Li2(z) = X^fc^i ^'^/^^ denotes a polylogarithm of 
second order. We now find for the imaginary part of the 
Chiral Magnetic conductivity at nonzero frequency and 
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momentum 



2 9 

e p--Po 



(q,+iM.)/T' 



^ 167r 

+ T ^si|2rT Lis 

s.i ,r— lb 

+ p [sgn(po)^(Po - P') + 7- - Po)] X 
log 1 + e 



(9,+tM.)/T 



}}. (43) 



In the limit of homogeneous magnetic fields (p ^ 0), 
the imaginary part of the Chiral Magnetic conductivity 
becomes 



967r 



s,t=± 



dg 



n{q + tfXs 



(44) 



g=|cj|/2 



The derivative with respect to the Fermi-Dirac distribu- 
tion shows that for small temperatures only states near 
the Fermi-surface contribute to the Chiral Magnetic con- 
ductivity. 

At zero temperature the imaginary part of the leading 
order Chiral Magnetic conductivity becomes 



967r 



st6{uj/2 + tfi,). (45) 

s,t=± 



The last equation shows that in the limit of p — s- the 
leading order contribution to the Chiral Magnetic con- 
ductivity develops resonances at w = ±2/1^. As argued in 
the previous section, these resonances can be attributed 
to particle-antiparticle pair production in the time de- 
pendent electromagnetic field. 

For large temperatures {T > fi^) we can approximate 
Eq. (m by 



g|u;|/T_g|^|/(2T) 



/is. (46) 



The imaginary part at large temperatures has a maxi- 
mum at w/T « 5.406 with value a" « 0.394cro- 



C. Real part 

At zero temperature we can now recover the real part 
of the Chiral Magnetic conductivity by applying the 
Kramers-Kronig relation, Eq. (|20p . We find 



<(^-0) 
<(^y^0) 



SfJ-s 



(47) 
(48) 



g^x(^) 
Co 





\ 1 1 

i < 












1 < 


1 


\ - 

1 1 



1 2 3 4 5 

FIG. 1: Real (red, solid) and imaginary (blue, dashed) part 
of the leading order Chiral Magnetic conductivity as a func- 
tion of frequency, at T = and p = for p = O.l/is. 
The result is normalized to the zero frequency conductivity 
(70 = eV5/(27r2). 



g^x(^) 



— ' j\ 












< \ 





















2 4 6 8 

FIG. 2: Real (red, solid) and imaginary (blue, dashed) part of 
the leading order normalized Chiral Magnetic conductivity as 
a function of frequency, ai T = 0, /j, — l.Sps and p = O.l/is. 
The result is normalized to the zero frequency conductivity 
cro = eV5/(27r2). 



For finite temperatures we were unfortunately not able 
to obtain an analytic expression for the real part. 

The delta function of the imaginary part at w = 
only contributes to 2/3 of the real part at w = 0. The 
other 1/3 part comes from the pair-production processes. 
This conclusion holds for any temperature. This is be- 
cause the both the real (see Eq. (501)) and imaginary 
part (Eq. (|ii)) ) at = are independent of tempera- 
ture. From the Kramers-Kronig relation it follows that 
for any nonzero frequency the imaginary part at w = 
does not contribute to the real part at a; 7^ 0. Therefore 
for any temperature the real part of the Chiral Magnetic 
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FIG. 3: Real (red, solid) and imaginary (blue, dashed) part of 
the leading order normalized Chiral Magnetic conductivity at 
high temperatures (T > fi^) for homogeneous magnetic fields 
(p = 0). At a; = the normalized conductivity is equal to 1. 



FIG. 4: Induced current in time-dependent magnetic field, 
Eq. (|49|l . as a function of time, at very high temperature. 
The results are plotted for different values of the characteristic 
time scale r of the magnetic field. 



conductivity drops from co at lu 
from LU — 0. 



to (7o/3 just away 



D. Discussion 

We display the real and imaginary part for T — 0, p — 
0.1^5 and ^ = in Fig. [1] As was argued at the end of 
the previous subsection, it can be seen in this figure that 
the real part of the Chiral Magnetic conductivity drops 
from (To at a; = to <7q/3 just away from lu ~ 0. Also 
the resonance at a; = 2/^5 is clearly visible. The width of 
the imaginary part at the resonance is equal to 2p. The 
real part of the conductivity becomes negative above the 
resonance frequency. This is a typical resonance behavior 
and implies that when the imaginary part vanishes the 
response is 180 degrees out of phase with the applied 
magnetic field. 

In Fig. [5] we display the real and imaginary part for 
T = 0, p — 0.1/i5 and fi — l.bfi^. In this case there are 
resonances at w = 5/X5 and w = /X5. Equation (|45)) shows 
that the imaginary part is proportional to cj^, therefore 
the second resonance at lu — 6/^5 is much stronger than 
the first one at tj = /is. Because the second resonance 
is due to the right-handed modes, and the first one due 
to left-handed, the contribution of the second resonance 
has opposite sign to the first resonance. 

The real and imaginary part of the Chiral Magnetic 
conductivity at high temperatures (T > /is) are displayed 
in Fig. [31 This figure is the most relevant for QCD at 
very high temperatures, since then loop corrections will 
be small. As argued in the previous subsection it can be 
seen in the figure that the real part of the conductivity 
drops from ctq at w = to cto/S just away from tu = 0. 

Let us now study the induced current in a magnetic 



field of the form created during heavy ion collisions. For 
simplicity we approximate the two colliding nuclei by 
point like particles like in Ref. 19|]. This gives a reason- 
able approximation to the more accurate methods dis- 
cussed in Refs. [13, [S] and is most reliable for large im- 
pact parameters. The magnetic field at the center of the 
collision can then be written as 



Bit) = 



1 



Bo, 



(49) 



with T = 6/(2sinhF) and cBq = SZ^em sinhy/6^. Here 
b denotes the impact parameter, Z the charge of the nu- 
cleus, and Y the beam rapidity. For Gold-Gold {Z — 79) 
collisions at 100 GeV per nucleon one has Y — 5.36. At 
typical large impact parameters (say 6 = 10 fm) one finds 
eBo ~ 1.9 X 10^ MeV^ and r = 0.05 fm/c. For 31 GeV 
per nucleon (Y = 4.19) Gold-Gold collisions one finds at 
& = 10 fm, eBo ~ 5.9 x lO" MeV^ and r = 0.15 fm/c. 
The Fourier transform of Eq. equals 



(50) 



B{lu) = 2t^\uj\Ki{t\uj\)Bo, 



where Ki{z) denotes the first-order modified Bessel func- 
tion of the second kind. 

For illustration purposes we will assume that our mag- 
netic field is (unlike in heavy ion collisions) homogeneous. 
The induced current can be found by applying Eq. (I19|) . 
We display the induced current in the magnetic field of 
Eq. (|49p in a system with nonzero chirality at very high 
temperatures in Fig. 21 The induced current is plotted as 
a function of time for three different characteristic time 
scales T of the magnetic field. 

In any general decaying magnetic field, the only rel- 
evant frequencies are the ones which are smaller than 
of order the inverse life-time of the magnetic field, lu ^ 
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0{1/t). In Fig. [3] it can be seen that even in the non- 
interacting case there still is sizable response at high tem- 
peratures as long as T > 1/(5T). Hence even for such fast 
changing fields there will really be a induced current, al- 
though it will be about 1/3 of the adiabatic approxima- 
tion which is j{t) — aoB(t). This can be more clearly 
seen in Fig. [H The red solid line describes a relatively 
slowly varying field with r = l/T, in that case indeed 
j{t) ~ aoB{t)/3. The red solid line can therefore also 
be seen as a guide to eye of the time dependence of the 
magnetic field. Realizing this, in a fast changing mag- 
netic field with r = O.l/T it clearly takes some time for 
the current to respond. The maximal current will be 
smaller, but it takes also more time for the current to 
diminish after the magnetic field has gone away. At late 
times the current can even become negative. The situ- 
ation in which r = O.l/T approximately corresponds to 
100 GeV per nuclcon Gold-Gold collisions at & = 10 fm. 

Thermal fluctuations will increase in magnitude when 
the temperature is increased. These fluctuations can 
cause the spins of the particles to align along the mag- 
netic field. Hence one would expect that when keeping r 
fixed and increasing the temperature the system will re- 
spond faster to the changing magnetic field. This can be 
seen in Fig. [U by increasing the temperature at fixed r 
one goes from the dotted to the solid line, which indicates 
faster response. 

The discussion will alter when loop corrections due 
to interactions are taken into account. In that case the 
quasi-particles obtain a thermal width F which is of order 
aT. Because of the thermal broadening we expect that 
the peak in the Chiral Magnetic conductivity (Fig. [3]) at 
uj = will get a width of order F. At the same time 
the value at w = will not change since it is constrained 
by the anomaly. As long as r ^ ©(l/F) the zero fre- 
quency result will therefore be a reasonable estimate and 
the induced current will be more or less equal to the adi- 
abatic approximation j{t) ^ aoB(t). Hence the stronger 
the system interacts at fixed fi^ and T, the stronger the 
response at small frequencies. Qualitatively one would 
expect such behavior, because the stronger a particle in- 
teracts the faster the particle aligns its spin and hence 
momentum along the field. 

Let us finally discuss how much charge is separated by 
the current. Let us divide our infinite system into an 
upper half (u) and a lower half (?), and set the boundary 
to be the plane z = 0. Using the the fact that the current 
is conserved = one finds that the change in charge 
in the upper hemisphere per unit of time equals 

^ = -|^d3xV.j= ld'xMz = 0,t). (51) 

By integrating this equation over time and using Eq. (|19|) 
we find that the change in charge in the upper hemisphere 
AQu equals 

/oo 
dtj{t)^aoAB{uj^O), (52) 
-oo 



where A is the area of the plane z = 0. Because of 
global charge conservation the charge change in the lower 
hemisphere is equal to AQi — —AQu- In Eq. ([52|) we 
find the surprising result that as long as linear response 
is valid and the chirality is constant the total induced 
charge difference between upper and lower hemisphere is 
independent of detailed dynamics and just determined 
by the zero frequency chiral magnetic conductivity (Tq. 
In the magnetic field of Eq. (|49p we obtain for the in- 
duced charge difference between the upper and lower 
hemisphere 

AQ EE AQu - AQi ^ AaoArBo. (53) 



VI. CONCLUSIONS 

In a quark gluon plasma it is possible to generate 
nonzero chirality by gluon configurations with nonzero 
topological charge. In the presence of a magnetic field 
nonzero chirality leads to a current along the field. This 
is the Chiral Magnetic Effect which can potentially give 
rise to observable effects in heavy ion collisions. Since the 
magnetic field in heavy ion collisions is rapidly decreasing 
as a function of time, it is desirable to study the Chiral 
Magnetic Effect in a time-dependent magnetic field. In 
this article, we have shown for the first time that such 
study is possible in a systematic way using linear response 
theory. 

To obtain the induced current in a time-dependent 
magnetic field we have derived a general Kubo-formula 
for the Chiral Magnetic conductivity. We have shown 
that the Chiral Magnetic conductivity is proportional to 
the antisymmetric part of the off-diagonal photon polar- 
ization tensor. 

Since we have applied linear response theory, our re- 
sults are only valid for small magnetic fields. This means 
that the magnetic fields should not alter the plasma dy- 
namics much, implying that the magnetic length l/y/eB 
has to be larger than the (color) electric screening length 
~ l/gT. Hence for reliable results the magnetic field 
should satisfy eB < g^T"^. 

We have computed the leading order Chiral Magnetic 
conductivity for constant chirality using perturbation 
theory. As such our result is only applicable for QCD 
at high temperatures where loop corrections can be ne- 
glected. We have shown that pair production in the time- 
dependent electromagnetic field gives rise to non-trivial 
behavior of the Chiral Magnetic conductivity at nonzero 
frequencies. Our result can be systematically improved 
by including loop corrections, in a similar way to what 
has been done for the electrical conductivity in for exam- 
ple Ref. [1^. 

The general formula for the Chiral Magnetic conduc- 
tivity we have obtained allows for the evaluation using 
other methods. For example as suggested in Ref. 40] one 
could study the Chiral Magnetic Effect in holographic 
models of QCD. Also since ^5 does not give rise to a 
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sign problem [l^l, one could in principle study the Chi- 
ral Magnetic conductivity using lattice QCD (see also 
Ref. [23] for another approach). However, because one 
has to perform an analytic continuation in order to obtain 
the retarded correlator this is not completely straightfor- 
ward to do. 

The main message of this calculation is the finding that 
even in the leading order the Chiral Magnetic conductiv- 
ity has sizable response at nonzero frequencies. This is 
a clear proof that even a rapidly decaying magnetic can 
give rise to a non- negligible current. By taking into ac- 
count interactions the quasi-particles will obtain a ther- 
mal width, which as we argued will increase the response 
at small frequencies but does not change the zero fre- 
quency result. 

This calculation of the Chiral Magnetic conductivity 
is a small step in order to improve our understanding 
of the dynamics of the Chiral Magnetic Effect. To ap- 
ply our results to heavy ion phenomenology one has to 
take into account also other dynamical effects like the 
time-dependence of the chirality, the radial flow, and 



possible screening mechanisms. Also it would be inter- 
esting to study the effects of a time-dependent magnetic 
field on other physical quantities and effects, like for ex- 
ample the chiral condensate, the chiral phase transition 
and dynamical chiral symmetry breaking. So far these 
have only been investigated in a constant magnetic field 

[n, El Si, s, Em, 113. 
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